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l/^ ' In this talk I discuss three main topics concerning the theoretical de- 

scription and observable signatures of possible phase transitions in nuclear 
collisions. The first one is related to the multifragmentation of equilibrated 
' sources and its connection to a liquid-gas phase transition in finite systems. 

*^ ■ The second one is dealing with the Coulomb excitation of ultrarelativistic 

, heavy ions resulting in their deep disintegration. The third topic is devoted 

to the description of a first order phase transition in rapidly expanding mat- 
ter. The resulting picture is that a strong collective fiow of matter will lead 
, to the fragmentation of a metastable phase into droplets. If the transition 

from quark- gluon plasma to hadron gas is of the first order, it will manifest 
itself by strong nonstatistical fiuctuations in observable hadron distributions. 

^ INTRODUCTION 

^ i A general goal of present and future experiments with heavy-ion beams is to 

I study the properties of strongly interacting matter away from the nuclear ground 

state. The main interest is focussed on searching for and studying possible phase 
transitions. Several phase transitions are predicted in different domains of tem- 
perature T and baryon density pb- There is no doubt that there should be a first 
order phase transition of the liquid-gas type in normal nuclear matter. This fol- 
lows simply from the existence of the nuclear bound state at the saturation density 
Po ~ 0.15 fm~^. Therefore, at ps < Po and low temperatures, T < Tc ~ 10 MeV, 
the matter will organize itself in the form of a mixed phase with droplets of nuclear 
liquid surrounded by the nucleon gas. The only problem is whether relatively small 
amounts of excited nuclear matter produced in nuclear collisions and its limited 
lifetime are sufficient to observe this phase transition. Based on recent data on the 
nuclear caloric curve |jl| and temperature fluctuations [0] I am tempting to give a 
positive answer to this question. This topic will be discussed in the first part of the 
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talk after a short description of the Statistical Multifragmentation Model (SMM) 
[0, |] which provides a basis for theoretical analysis. 

The situation at high T and nonzero baryon chemical potential fi {ps > 0) is not 
so clear, although everybody is sure that the deconfinement and chiral transitions 
should occur somewhere. The phase structure of QCD is not yet fully understood. 
Reliable lattice calculations exist only for n = (pb = 0) where they predict a 
second order phase transition or crossover at T ^ 160 MeV. As model calculations 
show, the phase diagram in the (T, p) plane may contain a first order transition line 
(below called the critical line) terminated at a (tri)critical point Possible 
signatures of this point in heavy-ion collisions are discussed in ref. 0. Under 
certain non-equilibrium conditions, a first order transition is also predicted for 
symmetric quark-antiquark matter f^. 

A striking feature of central heavy-ion collisions at high energies, confirmed in 
many experiments (see e.g. 0, |lOl), is a very strong collective expansion of matter. 
The applicability of equilibrium concepts for describing phase transitions under 
such conditions becomes questionable. In the last part of the talk I demonstrate 
that non-equilibrium phase transitions in rapidly expanding matter can lead to 



interesting phenomena which, in a certain sense, can be even easier to observe [11 



In the middle part of the talk I address a question which is closely related to 
the main topic of this conference. It illustrates how the knowledge accumulated in 
intermediate-energy heavy-ion physics can be used for ultrarelativistic heavy-ion 
colliders. Namely, I will discuss the excitation of nuclei by Lorentz-contracted and 
strongly-enhanced Coulomb fields of ultrarelativistic heavy ions. As well known, 
this process can be treated in terms of equivalent photons. Their flux grows linearly 
with the squared nuclear charge, and their characteristic energy is proportional to 
the relative Lorentz factor of colliding nuclei. This is why this Coulomb excitation 
of nuclei becomes especially important in high-energy heavy-ion colliders such as 
RHIC and LHC. The calculations P^ , [T^ show that in such colliders the equiva- 
lent photon spectrum extends far above the Giant Resonance region, into the GeV 
domain. The absorption of such a photon by a nucleus leads to its high excitation 
and subsequent disintegration. This might be an important factor determining a 
lifetime of ultrarelativistic heavy-ion beams. 



STATISTICAL MULTIFRAGMENTATION AND 
LIQUID-GAS PHASE TRANSITION 

When a nucleus is suddenly heated up to a temperature T it starts expanding 
to adjust a new equilibrium density Po(^) which is less than the equilibrium density 
at zero temperature po- If the initial temperature is high enough the expansion 
is unlimited. At some stage of expansion the system enters into the spinodal 
region, where the homogeneous distribution of matter becomes thermodynamically 
unstable. Therefore, the nucleons form smaller and bigger clusters or droplets 
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with density close to po- This clusterization process resembles a hquid-gas phase 
transition in ordinary fluids. In the transition region the matter is very soft in a 
sense that the sound velocity is close to zero (soft point). This means also that 
the expansion is slow and the system has enough time to find a most favorable 
cluster-size distribution maximizing the entropy. 

At a later stage of expansion the system reaches a so-called freeze-out state 
when clusters cease to interact with each other. This break-up state of the system 
can be described within a statistical approach. In 1985 we have constructed a 
Statistical Multifragmentation Model (SMM) 0, ^ which up to now is one of the 
most successful realizations of this approach for finite nuclear systems. The model 
and its numerous applications are described in detail in a recent review [|14|. A 



similar model was also constructed by Gross [|T5l. In this talk I outline only some 



general features of the SMM and give a few examples of how it works. 

It is assumed that at break-up the system consists of primary hot fragments 
and nucleons in thermal equilibrium. Each break-up channel or partition, /, is 
specified by the multiplicities of different species, Nazi constrained by the total 
baryon number and charge Zq. The total fragment multiplicity is defined as 
M = J2az^az- The probabilities of different break-up channels are calculated in 
an approximate microcanonical way according to their statistical weights, 

Wf^exp[Sf{E*,V,Ao,Zo)], (1) 

where Sf is the entropy of a channel / at excitation energy E* and volume V. 

Translational degrees of freedom of fragments are described by the Boltzmann 
statistics while the internal excitations of individual fragments with A > 4 are cal- 
culated according to the quantum liquid-drop model. An ensemble of microscopic 
states corresponding to a break-up channel / is characterized by a temperature Tf 
which is determined from the energy balance equation 

|t(M - 1) + E Eaz{T)Naz + Ef{V) -Qf = E* . (2) 

Here the first term comes from the translational motion, the second term includes 
internal excitation energies of individual fragments, the third term is the Coulomb 
interaction energy and the last one is the Q-value of the channel /. The exci- 
tation energy E* is measured with respect to the ground state of the compound 
nucleus {Ao,Zq). It is fixed for all fragmentation channels while the temperature 
Tf fluctuates from channel to channel. 

The total break-up volume is parametrized as V = {1 + k)Vo, where Vq is the 
compound nucleus volume at normal density and the model parameter k is the 
same for all channels. The entropy associated with the translational motion of 
fragments is determined by the "free" volume, Vf, which is only a fraction of the 
total break-up volume V. In the SMM Vf{M) is parametrized in such a way that 
it grows almost linearly with the primary fragment multiplicity M or equivalently, 
with the excitation energy e* = E*/Aq of the system [jl4| . 
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At given inputs Aq, Zq and e* the individual multifragment configurations are 
generated by the Monte Carlo method. After the break-up the hot primary frag- 
ments propagate in a common Coulomb field and loose their excitation. The most 
important de-excitation mechanisms included in the SMM are the simultane- 
ous Fermi break-up of lighter fragments {A < 16) and the evaporation from heavier 
fragments, including the compound-like residues. In refs. |T7| , |TB|, [T^, ^ ^ 



one can find fresh examples showing how well the SMM works in describing the 
multifragmentation of thermalized sources. 



In ref. |T6[ an equation of state of a multifragment system was calculated for the 



grand canonical version of the SMM. As expected, it shows clear signs of a liquid- 
gas phase transition with a critical temperature of about 7 MeV. In the transition 
region the pressure isotherms are very flat indicating that the sound velocity is 
very small. In this region the compressibility and specific heat have nonmonotonic 
behaviour. 

The most interesting prediction of the statistical model, a plateau in the caloric 
curve T{e*), was formulated already in 1985 |2^. Since that time it was a challenge 
for experimentalists to measure the nuclear caloric curve. First measurements 
were performed at GSI by the ALADIN collaboration only in 1995 They 
have shown an impressive agreement with the theoretical prediction. These results 
have initiated an avalanche of other measurements and a lively discussion in the 
community (see latest ALADIN results in ref. p5|). 



Most temperature measurements are based on the Albergo method |26| relat- 
ing the temperature to the double ratio of isotope yields. The analysis shows (see 
for instance |1], |2^) that the temperatures extracted by this method are very sen- 
sitive to the side-feeding and nuclear structure effects. According to SMM the 
observed light isotopes are produced mainly by the secondary decays of hot pri- 
mary fragments. This leads to a difference between the isotopic temperatures and 
true thermodynamical temperature at freeze-out (see detailed analysis in ref. [|I|). 
In particular, isotopic temperatures have typically a less pronounced plateau than 
the true temperature, which can even have a backbending. In recent years several 
comparisons have been made (see examples in refs. 0, 0, ^ which generally show 
very good agreement between the theory and experiment. 

One should bear in mind that the ALADIN caloric curves are measured for a 
wide ensemble of decaying sources associated with the projectile or target specta- 
tors produced in peripheral nuclear collisions. The thermodynamical significance of 
such observations would increase if the measurements were done for a fixed source 
size with a varying excitation energy. Also the temperature measurements on the 
event-by-event basis would make it possible to study its fluctuations and therefore 
the heat capacity of the nuclear system. Such an analysis was performed recently 
by the Bologna group in the study of quasi-projectile (QP) fragmentation in 
peripheral Au+Au collisions at 35 A MeV. In this analysis only the events with 
reconstructed QP charges 70 < Zqp < 88 were included. The excitation energy, 
determined by a calorimetric method, varied for these events from 0.5 to about 8 



4 




Figure 1: Experimental (a) and SMM generated (c) Campi scatter plots Three cuts are 
introduced to select liquid-like events (Cut 1), gas-like events (Cut 3) and critical events 
(Cut 2). Panels b) and d) show the correlation between temperatures and excitation 
energies for experimental and SMM generated events. Events belonging to different cuts 
are shown by blue (1), red (2) and green (3) squares. Sizes of the squares are proportional 
to the yields. The solid lines in panels b) and d) show the mean temperatures for all the 
events at given e*. 
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Mev/nucleon. 

Measuring temperatures event by event is of course a nontrivial task. An at- 
tempt of estimating the event "temperature" (below denoted by 9) was made in ref. 
[0]. The idea is to apply the energy balance equation which is used in the SMM, 
but now for the experimental events. Of course, this requires certain assumptions 
on how the observed partitions, involving cold reaction products, are related to the 
original partitions consisting of hot primary fragments. Therefore, it was assumed 
that the light particles detected in a partition were produced by de-excitation of 
hot primary fragments. For reconstructing a primary partition these light par- 
ticles were shared among the detected fragments proportionally to their charges 
and assuming the charge-to-mass ratio as in the entrance channel. Applying this 
procedure for the asymptotic SMM events showed that the correlation between the 
microcanonical temperature and excitation energy was reproduced within 5%. 

Fig. 1 shows the scatter plots in the (T, e*) plane for experimental (b) and 
SMM generated (d) events. The ensemble-averaged temperatures are indicated by 
the solid lines. Their behaviour is typical for caloric curves measured by other 
methods. In addition to the flattening of the average temperature at T 6 MeV, 
one can clearly see the broadening of the distributions in the transition region at 

= 4 ^ 3 MeV/nucleon. The quantity characterizing energy fluctuations is heat 
capacity. For a canonical ensemble at constant volume it can be expressed as 

- 7^ - 7^ ■ {'i) 

It is not clear whether the constant volume condition applies to actual freeze-out 
configurations but nevertheless studying the energy fluctuations provides an addi- 
tional and important information compared to the average characteristics. Indeed, 
applying Eq. (||) for scatter plots of Fig. 1 reveals a peak in Cy at temperatures 
around 6 MeV. This behaviour was also predicted theoretically a long time ago 
[SHI- 

Another way of characterizing the critical behaviour is to analyze the conditional 
moments of fragment multiplicity distributions introduced by Campi [^. Fig. 1 



a) shows, for each event j, the experimental correlation between the logarithm of 

■'big J 



the charge of the largest fragment, ln(Z^^^), and the logarithm of the correspond- 



ing second moment of the multiplicity distribution. In {1112'^) (Campi scatter plot). 
Fig. 1 c) shows the same for events generated by the SMM. As expected for a 
system experiencing a phase transition these plots exhibit two branches: an upper 
branch with an average negative slop, corresponding to under-critical events, and 
a lower branch with a positive slop that corresponds to super-critical events. The 
two branches meet in a central region signalling the approach to a critical point. 
This trend is nicely reproduced in Fig. 1 by both the experiment and the theory. 

We have made three cuts in these scatter plots selecting the upper branch (Cut 
1), the lower branch (Cut 3) and the central region (Cut 2) and analyzed the 
events falling in each of the three zones. The fragment charge distributions in 
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these three zones exhibit shapes going from a U-shape in Cut 1, characteristic of 
the evaporation events at low excitation energies, to an exponential one in Cut 3, 
characteristic of the vaporization events at high excitations. In Cut 2 a power-low 
fragment charge distribution Z""^ with r ^ 2.2 is observed as expected according 
to the Fisher's droplet model for fragment formation near the critical point of a 
liquid-gas phase transition (see also an interesting analysis of ref. [P^]). 

The contributions of these three types of events to the caloric curves are shown 
in Fig. 1 for experiment (b) and for theory (d). It is clearly seen for both the 
data and the SMM, that in Cuts 1 and 3 besides normal events there are un- 
usual events (although with low probability) which lie far from the average T{e*) 
behaviour. These are compound-like states with very high temperatures and vapor- 
ization events with low temperatures. For these events one can make the analogy 
respectively with an overheated liquid and a super-cooled gas in the ordinary liquid- 
gas phase transition. Here we see the advantage of a finite system where not only 
the most probable states but also the metastable states can be produced with a 
finite probability. In my opinion, the observation of these metastable states is the 
best indication that we are dealing here with the first order phase transition of the 



liquid-gas type. These interesting questions were further studied in ref. [31 



ELECTROMAGNETIC EXCITATION OF 
ULTRARELATIVISTIC HEAVY IONS 



It has become clear in recent years [12, M that high nuclear excitations can 



be induced by the Coulomb fields of ultrarelativistic heavy ions. Following the 
famous Waizsacker- Williams method the Lorentz contracted Coulomb field of an 
ultrarelativistic projectile in the rest frame of a target nucleus (and vice versa) can 
be represented as a beam of equivalent or virtual photons. The flux of equivalent 
photons with energy in a collision of nuclei with charge Z at impact parameter 
b is given by the standard formula 



N{E^,b) 



7r2 (3^E^b^ 



(4) 



where a is the fine structure constant, jS = v/c and 7 = a/1 — j3'^ is the relative 
Lorentz factor. The variable x in the modified Bessel functions Kq^i{x) is defined 
as a; = E^b / {P'yhc) . Since Kq i drop exponentially at large arguments, the main 
contribution to the virtual photon flux comes from the region a; ~ 1. Thus the 
characteristic energy of virtual photons grows linearly with 7. This explains why 
the relativistic Coulomb excitation is very important for ultrarelativistic heavy-ion 
beams where both 7 and Z are large. For colliding beams 7 = 27^^^^ — 1 that 
gives 2 ■ 10^ and 10'' for RHIC and LHC respectively. This brings the spectrum of 
virtual photons into the GeV energy domain, i. e. much above the traditionally 
studied Giant Resonance (GR) and Delta-resonance regions. The absorption of 
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Figure 2: Theoretical predictions for inclusive cross sections for emitting nucleons, 
pions, intermediate mass fragments (IMF: 3 < Z < 30) and fission fragments (FIS: 
30 < Z < 50) in the electromagnetic dissociation of Au nuclei as functions of the cm. 
energy ^/s (the SPS, RHIC and LHC energies are indicated by arrows). The thick dotted 
line shows the total ED cross section for Au beams. 



such a high-energy photon will result in a very high nuclear excitation sufficient 
for its total disintegration. 

In ref. the description of nuclear photoabsorption was extended to the 



photon energies much above the GR region, where the excitation of individual nu- 
cleons and multiple pion production are the dominant reaction channels. A model 
of electromagnetic dissociation (ED) taking into account these high-energy photon 
absorption channels was constructed in ref. |jl2|]. According to this model, the 
fast hadrons produced after the photon absorption initiate a cascade of subsequent 
collisions with the intranuclear nucleons leading to the fast particle emission and 
heating of a residual nucleus. This stage is described by the Intranuclear Cascade 
Model (INC). At a later stage the nucleus undergoes de-excitation by means of 
the evaporation of nucleons and lightest fragments, binary fission or multifragmen- 
tation. The latter process becomes important at ultrarelativistic beam energies, 
when the excitation energy of residual nuclei exceeds 3-4 MeV/nucleon. This stage 
of the reaction is described by the SMM. 
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To include all the processes described above, in ref. we have developed a 
specialized computer code RELDIS aimed at the Monte Carlo simulation of the 
Relativistic ELectromagnetic Dissociation of nuclei. The simulation begins with 
generating the single- or double-photon absorption process. Then the INC model is 
used to calculate the fast particle emission and the characteristics of residual nuclei. 
Finally, de-excitation of thermalized residual nuclei is simulated by the SMM. 

The cross section of the photo-nuclear (7A) reaction induced by a photon of 
energy is expressed as 



dE^ Jb. 



00 



"^""^^ cT^AiE^) I N{E^,b)2TTbdb, (5) 



where bmm ~ {Rp + Rt) is a minimal impact parameter for heavy-ion collisions 
without nuclear overlap, a^A^E^) is an appropriate photo-absorption cross section, 
either measured for the A-nucleus with real photons or calculated within a model. 
The total ED cross section, aED, is obtained by integrating Eq. (|) by dE^ from to 



00. The calculations [|12| show that the total ED cross sections for RHIC and LHC 
are very large, 100 b and 200 b respectively. Accordingly, the ED reaction rates are 
much higher than those for nuclear interactions, although the ED events are much 
less violent. For instance, at expected RHIC luminosity L ^ 10^^ cm^^s^^ the ED 
reaction rate will be 10^ interactions per second. Together with the electron capture 
reactions the ED processes will be the important factors reducing the lifetime of 
ultrarelativistic heavy-ion beams compared with the proton ones. 

We have applied the RELDIS code for calculating the ED characteristics for 
several heavy-ion beams. The model is in a reasonable agreement with experimental 
data, when available. We have also made predictions for the reactions: 160A 
GeV Pb+Pb (SPS), lOOA+lOOA GeV Au+Au (RHIC) and 2.75A+2.75A TeV 
Pb+Pb (LHC). The inclusive (multiplicity weighted) cross sections for emitting 
nucleons, pions and nuclear fragments in the electromagnetic dissociation of one 
of the colliding Au nuclei are shown in Fig. 2 as functions of the incident cm. 
energy. Nuclear fragments are divided in two groups: fission fragments (30<Z< 
50) and Intermediate Mass Fragments (IMFs, 3<Z<30), which are associated with 
the multifragmentation. One can clearly see a steep rise in the yields of all species, 
especially IMFs, when the incident energy grows from the SPS to RHIC and LHC 
domain. The inclusive cross section for neutron emission is especially large, above 
1000 b at RHIC and LHC. The average neutron multiplicities are predicted to be 
4.1, 7.2 and 8.8 at SPS, RHIC and LHC respectively. 

The predicted neutron multiplicity distributions are shown in Fig. 3. They have 
a nontrivial structure. There is a strong peak at In emission channel associated 
with the GR decay. On the other hand, there is a long tail of multiple neutron emis- 
sion associated with more violent reaction channels, from the direct knock-out and 
evaporation from the compound nucleus to fission and multifragmentation. This 
is where our model including all these channels shows its strength. One can see, 
for example, that the probability to emit more than 20 neutrons is quite noticeable 
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Figure 3: Normalized neutron multiplicity distributions for the electromagnetic dissoci- 
ation of Pb nuclei at LHC and SPS (solid and dotted histograms, respectively) and Au 
nuclei at RHIC energies (dashed histogram). Calculations are made with RELDIS code. 



(~ 5% at RHIC). These results might be important for designing neutron-sensitive 
zero-degree calorimeters at RHIC and LHC. One of such proposals was made re- 
cently in ref. but only the In channel was considered there. 



FIRST ORDER PHASE TRANSITION IN FAST 

DYNAMICS 

The implications of a strong collective expansion on the liquid-gas phase tran- 
sition were discussed in ref. [|^. Here I will focus on consequences of the strong 
collective flow of matter for a possible first order chiral transition. I will assume 
that the collective velocity field is described locally by the Hubble law, t>(r) = H -r, 
where the Hubble "constant" H may in general depends on time. 

To make the discussion more concrete, I adopt a picture of the chiral phase 
transition predicted by the linear sigma-model with constituent quarks pS. Then 



the mean chiral field $ = (a, vr) serves as an order parameter. The model respects 
chiral symmetry, which is spontaneously broken in the vacuum where a = fj^, 
71 = 0. The effective thermodynamic potential f2(T, //;$) depends, besides $, 
on temperature T and baryon chemical potential fi. The schematic behaviour of 
f2(T, n; $) as a function of the order parameter field a at vr = is shown in Fig. 4. 
The minima of Q correspond to the stable or metastable states of matter under the 
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Figure 4: Schematic view of the effective thermodynamic potential per volume Vt/V as 
a function of the order parameter field cj at tt = 0, as predicted by the linear cr- model in 
the chiral limit m.,^ = [35]. The curves from bottom to top correspond to the different 
stages of the isentropic expansion of homogeneous matter starting from r=100 MeV and 
//=750 MeV (curve 1). The upper curve 5 is the vacuum potential. The other curves are 
discussed in the text. 

condition of thermodynamical equilibrium, where the pressure is P = — fimin/^- 
The curves from bottom to top correspond to different stages of the isentropic 
expansion of homogeneous matter. Each curve represents a certain point on the 
(T, /i) trajectory. As one can see from the figure, the model of ref. reveals a 
rather weak first order phase transition, although some other models 0, ^ predict 
a stronger transition. The discussion below is quite general. 

Assume that at some early stage of the reaction the thermal equilibrium is es- 
tablished, and partonic matter is in a "high energy density" phase Q. This state 
corresponds to the absolute minimum of Vl with the order parameter close to zero, 
cr 0, vr ^ 0, and chiral symmetry restored (curve 1). Due to a very high internal 
pressure, Q matter will expand and cool down. At some stage a metastable mini- 
mum appears in at a finite value of a corresponding to a "low energy density" 
phase H, in which chiral symmetry is spontaneously broken. At some later time, 
the critical line in the (T, /i) plane is crossed where the Q and H minima have 
equal depths, i.e. Pr = Pq (curve 2). At later times the H phase becomes more 
favorable (curve 3), but the two phases are still separated by a potential barrier. 
If the expansion of the Q phase continues until the barrier vanishes (curve 4), the 
system will find itself in an absolutely unstable state at a maximum of the thermo- 
dynamic potential. Therefore, it will freely roll down into the lower energy state 
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corresponding to the H phase. This situation is known as a spinodal instabihty. 
As well known, a first order phase transition proceeds through the nucleation 



process. According to the standard theory of homogeneous nucleation super 



critical bubbles of the H phase appear only below the critical line, when > Pq- 
In rapidly expanding matter the nucleation picture might be very different. As 
shown in ref. |]TT|, the phase separation in this case can start as early as the 
metastable H state appears in the thermodynamic potential, and a stable interface 
between the two phases may exist. An appreciable amount of nucleation bubbles 
and even empty cavities may be created already above the critical line. 

The bubble formation and growth will also continue below the critical line. 
Previously formed bubbles will now grow faster due to increasing pressure differ- 
ence, Ph ~ -Pq > 0, between the two phases. It is most likely that the conversion 
of Q matter on the bubble boundary is not fast enough to saturate the H phase. 
Therefore, a fast expansion may lead to a deeper cooling of the H phase inside the 
bubbles compared to the surrounding Q matter. Strictly speaking, such a system 
cannot be characterized by the unique temperature. At some stage the H bubbles 
will percolate, and the topology of the system will change to isolated regions of the 
Q phase (Q droplets) surrounded by the undersaturated vapor of the H phase. 

The characteristic droplet size can be estimated by applying the energy balance 
consideration, proposed by Grady [^, ^ in the study of dynamical fragmentation 
of fluids. The idea is that the fragmentation of expanding matter is a local pro- 
cess minimizing the sum of surface and kinetic (dilational) energies per fragment 



volume. As shown in ref. this prescription works fairly well also for multifrag- 
mentation of expanding nuclei, where the standard statistical approach fails. 

Let us imagine an expanding spherical Q droplet of radius R, embedded in the 
background of the dilute H phase. The change of the thermodynamic potential, 
AQ, compared to the uniform H phase can be easily estimated within the thin- 
wall approximation |TT|. According to the Grady's prescription, the quantity to be 
minimized is AQ per droplet volume, V oc R^, that is 

(^) =-iP^-Pn) + 'i + ^ASH^R' . (6) 

V V /droplet rC iU 

Here AS = Sq—Syl is the difference of the bulk energy densities of the two phases, 7 
is the interface energy per unit area. One should notice that the last term, i.e. the 
change in the collective kinetic energy, is positive because £q> Sn. This term acts 
here as an effective long-range potential, similar to the Coulomb potential in nuclei. 
Since the bulk term does not depend on R the minimization condition constitutes 
the balance between the collective kinetic energy and the interface energy. This 
leads to an optimum droplet radius 

1/3 



R* 



(^) . (T) 



One can say that the metastable Q phase is torn apart by a mechanical strain 
associated with the collective expansion. This phenomenon has a direct analogy 
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with the fragmentation of pressurized fluids leaving nozzles PP| , In a similar 
way, splashed water forms droplets which have little to do with the equilibrium 
liquid-gas phase transition. 

In the lowest-order approximation the characteristic droplet mass can be cal- 
culated as M* A£V . It is natural to think that nucleons and heavy mesons are 
smallest droplets of the Q phase. For numerical estimates I take 7 = 10 MeV/fm^ 
and = 0.5 GeV/fm^, i.e. the energy density inside the nucleon. For the Hubble 
constant I consider two possibilities: H^^ = 20 fm/, representing a slow expansion 
from a soft point, and = 6 fm/c typical for a fast expansion. Substituting these 
values in Eq. (0) one gets i?*=3.4 fm and 1.5 fm for the slow and fast expansion 
respectively. These two values of R* give M* of about 100 GeV and 10 GeV, re- 
spectively. At ultrarelativistic energies the collective expansion is very anisotropic, 
with the strongest component along the beam axes. For the predominantly 1-d 
expansion one should expect the formation of slab-like structures with intermittent 
layers of Q and H phases. 

After separation the droplets recede from each other according to the global 
Hubble expansion, predominantly along the beam direction. Hence their center- 
of-mass rapidities are in one-to-one correspondence with their spatial positions. 
Presumably they will be distributed more or less evenly between the target and pro- 
jectile rapidities. Since rescatterings in the dilute H phase are rare, most hadrons 
produced from individual droplets will go directly into detectors. One can guess 
that the number of produced hadrons is proportional to the droplet mass. Each 
droplet will give a bump in the hadron rapidity distribution around its center-of- 
mass rapidity. If emitted particles have a Boltzmann spectrum, the width of the 



bump will be 6y ~ 2-y/T/m, where T is the droplet temperature and m is the par- 
ticle mass. At T ~ 100 MeV this gives 6y ^ 2 for pions and 6y ^ 1 for nucleons. 
These spectra might be slightly modified by the residual expansion of droplets and 
their transverse motion. The resulting rapidity distribution in a single event will be 
a superposition of contributions from different droplets, and therefore it will exhibit 
strong non-statistical fluctuations. The fluctuations will be more pronounced if pri- 
mordial droplets are big, as expected in the vicinity of the soft point. If droplets 
as heavy as 100 GeV are formed, each of them will emit up to ~300 pions within a 
narrow rapidity interval, 5y ^ 1. Such bumps can be easily resolved and analyzed. 
The fluctuations will be less pronounced if many small droplets shine in the same 
rapidity interval. Critical fluctuations of similar nature were discussed in ref. 



Some unusual events produced by high-energy cosmic nuclei have been already 
seen by the JACEE collaboration Unfortunately, they are very few and it is 
difficult to draw definite conclusions by analyzing them. We should be prepared 
to see plenty of such events in the future RHIC and LHC experiments. It is clear 
that the nontrivial structure of the hadronic spectra will be washed out to a great 
extent when averaging over many events. Therefore, more sophisticated meth- 
ods of the event sample analysis should be used. The simplest one is to search 
for non-statistical fiuctuations in the hadron multiplicity distributions measured 
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in a fixed rapidity bin One can also study the correlation of multiplicities 

in neighbouring rapidity bins, bump-bump correlations etc. Such standard meth- 
ods as intermittency and commulant moments wavelet transforms 0|, HBT 
interferometry |^ can also be useful. All these studies should be done at dif- 
ferent collision energies to identify the phase transition threshold. The predicted 
dependence on the Hubble constant and the reaction geometry can be checked in 
collisions with different ion masses and impact parameters. 



CONCLUSIONS 

• The statistical approach (SMM) works well in situations when thermalized 
sources are well defined and no significant collective flow is present. 

• The quantitative agreement of SMM with recent data on the caloric curve and 
temperature fluctuations provides a strong indication on the nuclear liquid- 
gas phase transition. The nuclear heat capacity has a peak at T ^ 6 MeV. 

• A first order phase transition in rapidly expanding matter should proceed 
through the nonequilibrium stage when a metastable phase splits into droplets. 
The primordial droplets should be biggest in the vicinity of a soft point when 
the expansion is slowest. 

• Hadron emission from droplets of the quark-gluon plasma should lead to large 
nonstatistical fluctuations in their rapidity spectra and multiplicity distribu- 
tions. The hadron abundances may reflect directly the chemical composition 
in the plasma phase. 

• Electromagnetic excitation of nuclei in ultrarelativistic heavy-ion colliders is 
an important reaction mechanism leading to the deep nuclear disintegration. 
The multiple neutron emission associated with this process may be used for 
monitoring ultrarelativistic heavy-ion beams. 

• And finally, we should use the lessons of the hquid-gas phase transition for 
future studies of the deconfinement-hadronization and chiral phase transitions 
in relativistic heavy-ion collisions. 
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